Direct estimation of functionals of density operators by local operations and classical 

communication 
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We present a method of determining important properties of a shared bipartite quantum state, 
within the "distant labs" paradigm, using only local operations and classical communication 
(LOCC). We apply this procedure to spectrum estimation of shared states, and locally imple- 
mentable structural physical approximations to incompletely positive maps. This procedure can 
also be applied to the estimation of channel capacity and measures of entanglement. 

PACS numbers: 03.67.-a,03.67.Hk,03.67.Mn 



There are many scenarios in quantum information sci- 
ence where it is necessary to estimate certain properties 
of a quantum state g, such as its spectrum, purity or de- 
gree of entanglement. Moreover, such estimations are of- 
ten needed when g is a bipartite state gAB , shared by two 
distant parties, Alice and Bob, who can perform only lo- 
cal operations and communicate classically (LOCC). The 
desired properties can be then estimated either by resort- 
ing to quantum state tomography Tl or more directly, 
e.g. via estimating non-linear functionals of gAB- The 
second method has the natural advantage of being more 
efficient, since we compute directly the desired properties 
without estimating any superfluous parameters. In fact 
the direct estimation has been successfully applied to lo- 
cal spectrum estimation [3], entanglement detection 0,01 
and the evaluation of one-qubit quantum channel capac- 
ities 0. However, the LOCC version of these techniques 
was left as an open problem. In this paper we show that 
the two basic techniques, namely, the estimation of non- 
linear functionals of quantum states and constructions of 
Structural Physical Approximations 13, admit LOCC 
implementation. This opens the possibility of the direct 
estimation of entanglement and some channel capacities 
using only LOCC. 

As a general remark, let us recall that a quantum op- 
eration A can be implemented using LOCC if it can be 
written as a convex sum 



A 



Pk Ak 



(1) 



where Ak acts on the subsystem at Alice's location and 
Bk on the subsystem at Bob's location, and represent 
the respective probabilities. 

Let us start with the estimation of non-linear func- 
tionals of gAB using quantum interferometry. Consider 
a typical set-up for single qubit interferometry, conve- 
niently expressed in terms of quantum gates and net- 
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FIG. 1: A modified Mach-Zender interferometer with coupling 
to an ancilla by a controlled-?/ gate. The interference pattern 
is modified by the factor ve'" = Tr [Up\. 



works: Hadamard gate, phase-shift gate, Hadamard gate, 
and measurement in the computational basis {|0), |1)}. 
We modify the interferometer by inserting a controUed- 
U operation between the Hadamard gates, with its con- 
trol on the qubit and with U acting on a quantum state 
P (Fig-CJ- The controlled- [/ models the interaction be- 
tween the qubit and an auxiliary system (ancilla), ini- 
tially in the state p, and it leads to modification of the ob- 
served interference pattern, by the factor ve*" = Tr \U p] . 
The factor v is the new visibility and a is the shift 
of the interference fringes, known as the Pancharatnam 
phase The observed modification of the fringes gives 
us an estimate of the average value of unitary operator 
U in state p 0. 

Suppose now that p is the quantum state of two sepa- 
rable subsystems, p — gA® gB and that we choose U to 
be the swap operator V, defined such that V\(t)) a\'4') b ~ 
|V')a|</')s, I?/'). In this case, the modification of 

the interference pattern will be v = Tr [^(giA ® Pb)] = 
Tr[pAPB], or the overlap between the input states gA 
and gs. If the two inputs states are equal, gA = gB ^ g, 
we obtain an estimation of the purity, Tr [gi^] . The gen- 
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FIG. 2: Network for remote estimation of non-linear function- 
als of bipartite density operators. Since Trfy'^^'gi®*] is real, 
Alice and Bob can omit their respective phase shifters. 



eralization of the swap operation V to the shift opera- 

tion (y«i</.i)i</)2)...|0fe) - \^km)...\(t>k-i), 

i = 1, ...,k), and the choice of p = g^^ as the input state, 
allows us to estimate multi-copy observables, Tr[£i'''], of 
an unknown state g 0; S 3 ■ 

Let us now extend this method to the LOCC scenario 
by constructing two local networks, one for Alice and one 
for Bob, in such a way that the global network is similar 
to the network with the controUed-shift. Unfortunately, 
the global shift operation V^'^^ cannot be implemented 
using only LOCC, since it does not admit decomposi- 
tion (Q. Thus, we will implement it indirectly, using the 
global network shown in Fig. [21 Alice and Bob share a 
number of copies of the state gAB £ B{TL'^). They group 
them respectively into sets of k elements, and run the 
local interferometric network on their respective halves 
of the state pab = Qab- ea,ch run of the experiment, 
they record and communicate their result. 

The individual interference patterns Alice and Bob 
record will depend only on their respective reduced den- 
sity operators. Alice will observe the visibility ya — 
Tr[g\] and Bob will observe the visibility vb = Ti[gg]. 
However, if they compare their individual observations, 
they will be able to extract information about the global 
density operator gAB, e.g. about 



TrigU-^r gTs [V^a^ 



(2) 



This is because Alice and Bob can estimate the probabili- 
ties Pij that in the measurement Alice's interfering qubit 
is found in state \i)A and Bob's in state \ j)A for i, j = 0, 1. 
These probabilities can be conveniently expressed as 



rTr 



(k) 

B 



(3) 



hence the formula for the basic non-linear functional of 
gAB reads 



Tr[e^ 



AB\ 



Poi 



Pll 



(4) 



In fact, the expression above is the expectation value 
{(Jz ®az), measured on Alice's and Bob's qubits (the two 
qubits that undergo interference). Given that we are able 
to directly estimate Tr[£i^^] for any integer value of fc, 
we can estimate the spectrum of gAB without resorting 
to a full state tomography. 

We next show how to implement Structural Physical 
Approximations within the LOCC constraint. Struc- 
tural Physical Approximations (SPAs) were introduced 
recently as tools for determining relevant parameters of 
density operators (see 0, 1^ for more details) . Basically 
the SPA of a mathematical operation A, denoted as A, is 
a physical operation, a process that can be carried out in 
a laboratory, that emulates the character of A. More pre- 
cisely, suppose A : B{H'^) i-^ BiTL'^) is a trace preserving 
map which does not represent any physical process, for 
example, an anti-unitary operation such as transposition. 
Then a convex sum 



A = aP+ (1 -q;)A, 



(5) 



where T) is the depolarizing map which sends any den- 
sity operator into the maximally mixed state, represents 
a physical process, i.e. a completely positive map, as 
long as a is sufficiently large. On top of this T), with 
its trivial structure, does not mask the structure of A. 
The Structural Physical Approximation to A is obtained 
by selecting, in the expression above, the threshold value 
a — {(PX)/{d?X + 1), where —A is the lowest eigenvalue 
of (1(g) A) pj.'''' and Pf^ is a maximally entangled state of 
a, d X d system 14] . 

Please note that the physical implementation of SPAs 
is not a trivial problem as the formula (j^J , which explic- 
itly contains the physically impossible map A, is of little 
guidance here. Let us also mention in passing that if A is 
not trace preserving then A may be implementable but 
only in a probabilistic sense e.g. via a post-selection. 

There are many examples of mathematical operations 
which although important in the formulation of the phys- 
ical theory do not represent a physical process. For exam- 
ple, mathematical criteria for entanglement involve pos- 
itive but not completely positive maps ^ and as such 
they are not directly implementable in a laboratory — 
they tacitly assume that a precise description of a quan- 
tum state of a physical system is given and that such 
operations are performed on the mathematical descrip- 
tion of the state rather than the system itself. 

If A does not represent any physical process then its 
trivial extension to a bipartite case, I® A, does not repre- 
sent a physical process either. Still, its SPA, I ® A, does 
describe a physical operation, but can it be implemented 
with LOCC? 
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The positive answer is obtained by putting I ® A into 
the tensor product form Let us start by writing it 
as 



)A = 01?®©+ (1 - q)I® A 
1-a 



{l~a + l3)I( 



+ {a -13) 



-P 



A 



f3 



^a- f3 a- P 
(l-a + /3)I®A+(a-/3)e 



V 



)0, 



where 



A 

e 



1 



I- a 
a 

a — [3 



V 



i-l). 



(6) 

(7) 

(8) 



Equation (jSJ docs not represent a convex sum of physi- 
cally implementablc maps for any values of a and /3 but 
if we choose 



/3 > (1 - a)X(f 
a > 



(9) 
(10) 



where —A is the minimum eigenvalue of I ® A(P^), then 

indeed I (g) A is a physical operation in the LOCC form. 
Note, however, that the map O is not trace preserving 
and as such it can be implemented only with a certain 
probability of success. The minimal parameters a and /3 
that satisfy inequalities Eqs. and 1(1171) are 



Ad* 



/3 = 



Ad* + 1 ' 
Ad* + 1 ■ 



(11) 
(12) 



Hence, the SPA I ® A can be implemented, by Alice and 
Bob, using only only LOCC. 

The SPAs have been employed to test for quantum 
entanglement 2|. Recall that a necessary and sufficient 
condition for a bi-partite state qab to be separable is I (8> 
Mqab) > 0, for all positive maps A This condition, 
when considering the SPA I A on qab , is equivalent to 



1(g) A 



QAB > 



d^A + r 



(13) 



where —A is the minimal eigenvalue of the state 
[(I(g)I) ® (I® A)] (P^') 0. Thus, by estimating the 
spectrum (or the lowest eigenvalue) of the state 



1(g) A 



QAB, we can directly detect quantum entangle- 



ment. Moreover, we have already shown that both I (g) A 
and the spectrum estimation of qab, via non- linear func- 
tional, can be performed using only LOCC, hence a 
direct detection of quantum entanglement within the 
LOCC scenario is also possible. 



Let us now comment briefly on other potential applica- 
tions of the methods presented above. Let a completely 
positive map A : B{'H'^) B{H'^) represent a quantum 
channel shared by Alice and Bob. An estimation of the 
channel capacity may involve either a channel tomogra- 
phy or a direct estimation. In the case of tomography 
Alice prepares a maximally entangled pair of particles in 
state and sends one half of the pair to Bob. They 
now share the state 



PA = [I g) A] P_f 



(14) 



From the Jamiolkowski isomorphism this bi-partite 
state encodes all properties of the channel A, so state 
tomography on qa is effectively channel tomography on 
A. However, given a bi-partite state q\, Alice and Bob 
can also use the LOCC techniques to directly estimate its 
desired properties. For example, it has been shown that 
a single qubit channel A has non-zero channel capacity 
if and only if the maximal eigenvalue of qa is strictly 
greater than i (see Q for details). This can be estimated 
directly via the spectrum estimation, which in the case 
of two qubits requires only 2x4 — 3 = 5 measurements of 
the type Cz^cr^ as opposed to the 15 parameters required 
for the state estimation. 

For Bell diagonal states (i.e. two-qubit states, whose 
eigenvectors are all maximally entangled), the entangle- 
ment of formation (or negativity, see below) can be in- 
ferred from its spectrum 0. Thus, if Alice and Bob 
share a Bell entangled state, they can estimate the de- 
gree of entanglement of their state through spectrum 
estimation only. An important subclass of Bell diago- 
nal states are the maximally correlated states, rank two 
states equivalent (up to Ua ® Ub transformations) to 
mixtures of two pure states, 1-0+) — "^(|0)|0) + |1)|1)) 
and = 7f(|0)|0) - |1)|1)). The one-way distill- 

able entanglement can be calculated for such states as 
= log 2 — S{q), which is a function solely of the spec- 
trum. Thus, instead of estimating the seven parameters 
required to describe maximally correlated states, we need 
only estimate five parameters. 

The estimation of entanglement measures (see [i3 for 
review) is known only for special cases, such as the 
computable measure of entanglement [ll|, N{qab) = 
log 1 1 1 1 — log(^- I Ail). This measure is valid for any 
shared bipartite state, with a maximally mixed reduced 
density operator of at least one sub-system, and it is 
a function of the spectrum {A^} of the partially trans- 
posed matrix q^^ = I g) T[qab), where T is the (in- 
completely positive) transposition map. Thus, we can 
estimate J^{qab) using only LOCC, if we choose A = T 
and then estimate the spectrum of the resultant state. 

Given any quantum channel A, A/'(pa) is the upper 
bound for one-way channel capacity. We obtain, there- 
fore, a necessary condition for non-zero one-way capacity 
Q^: if A/'(pa) = 0, the two-way channel capacity must 
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vanish '1T\ (this can be easily seen using distillation and 
binding entanglement channel Hence, the positive 

partial transpose (PPT) test is a strong necessary test of 
quantum non-zero channel capacity. 

To conclude we have demonstrated that both direct 
spectrum estimations and the structural physical approx- 
imations can be implemented in the case of bi-partite 
states using only local operations and classical communi- 
cation. This leads to more direct, LOCC type, detections 
and estimations of quantum entanglement and of some 
properties of quantum channels. Direct estimations of 
specific properties have the natural advantage over the 
state tomography because they avoid estimating super- 
fluous parameters. Still, the exact comparison of the use 
of physical resources in tomography and direct estima- 
tions depends very much on the physical implementations 
of these techniques. Our objective here is to provide ad- 
ditional tools for quantum information processing rather 
than comparative studies of these tools. 
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